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The main model studied in this paper is a lattice of nearest neighbors coupled 
pendula. For certain localized coupling we prove existence of energy transfer and 
estimate its speed. 

1 The description of the motion. 

We consider a system of pendula with a nearest neighbors coupling: 

d 

Xi + sin. Xi = -e— - P{xi- !,Xi,Xi +1 ,e), i £ Z, (1) 
axi 

where the interaction potential /3 is localized and will be defined later. This system 
can be written in the Hamiltonian form with x = {xi}i£%, y = {yi}i£%, Xi and 
Ui £ R, with the the Hamiltonian 



H e (x,y) = -^- + ( cos X{ - 1) + e/3(xi- 1 ,Xi,x i+1 ,e) 



iez, 



y + V{xi) + ep{x i _ 1 ,x i , x i+1 ,e), 



(2) 



where V{x) = — cos a; — 1 is the pendulum potential. 



The system is near-integrable for small e, and most (in the sense of measure) 
of the systems's phase space is taken up by invariant KAM tori. In particular, for 
most initial data the energy of each pendulum will stay close to its initial value for 
all time. Nevertheless, we will show that this is not so for some motions, where the 
energy can slowly "seep" from one pendulum to another. We will in fact prove that 
for an arbitrarily small e and for any sequence of integers a = (..., cr_i, oo, <Ji, . . .) 
such that do = 0, \<jj — <?j+i\ = 1 for all j £ Z there exists a sequence of times 
(. . . , t-i, to, ti, . . .) (depending on e) such that at time tk the a^-th pendulum has 
most of the system's energy. In particular, one can make the energy wander along 
the chain of the pendula in any prescribed fashion, advancing to the right any 
number of steps, retreating to the left any number of steps, and so on. 
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Figure 1: The "running" and the near-heteroclinic motion are the building blocks 
of the dynamics. 

From now on we fix the energy of the system to be 10- Below we shall concentrate 
on the case of periodic collection of 4 pendula, i.e. of the index i (mod 4). The 
proof in the general periodic case i € Z/pZ is quite similar and necessary remarks 
are made along the proof. 
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Figure 2: A mechanical interpretation of @ with j3 — sin(xj_i — 2xj + a^+i). 

We note as a side remark that the space discretization of the sin-Gordon equation 
u u — u ss — sin u results in a system of pendula with elastic coupling [BWEj IKP1 
IWES] : 

P(xj-i,Xj,x j+1 ) = a(xj_i +x j+ i ~ 2xj); (3) 

this corresponds to an elastic torsional coupling between the neighbors. In partic- 
ular, as the angle — Xj oo we have /3 — > oo. By contrast, the coupling we 
consider in this paper can be interpreted as coming from a spring connecting points 
on a circle with angular coordinates shown in Figure [2] 

The coupling in our system is, however, localized, as described next. The class of 
coupling functions /3 for which our results hold is defined as follows: Let ij : R + — > R 
be a C°° bump function: r)(x) > for \x\ < 1 and rj(x) = for |a;| > 1. Exact form 

1 Energy corresponds to all pendula upside down and at rest. Indeed, the maximum of the 
potential energy V(x) = — cosx — 1 of an individual pendulum is and is achieved at x = it, an 
upside-down position. 
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of rj is not important, and in particular, no monotonicity properties are assumed; 
in particular, 77 is allowed to have many local maxima and minima, as long as the 
above conditions hold. We now dehne 

(3(x,e)=e r v O X ~ 2nnl ), x = ( Xl ,x 2 ,x 3 )eR 3 , (4) 

and fix r > 3 from now on. This is a C°°-smooth 27r-periodic function in each 
Xj, j — 1,2,3, or, equivalently, a function on 27rM 3 /Z 3 . Note that the C r -norm 
of e/3(-,e) tends to zero as e — > 0, while the norms of order r + 2 and higher are 
unbounded for e — > 0. We fix any finite r > 3 from now on. 

We will sometimes refer to connected components of support of /3 as lenses: in 
fact, they act by defocusing geodesies in the Jacobi metric, as explained later, as in 

En]. 



According to the main theorem, stated next, the energy 

Ej :=^ + V( Xj ) (5) 

at the jth site can pass from one site to another in an arbitrarily prescribed sequence 
of steps, as illustrated in Figure [3] Here is a more precise statement. 




• • • • • • » 

-3-2-1 1 2 3 



Figure 3: Any path in the graph Z can be shadowed by a solution of ((XJ) - 

Theorem 1. Let us fix the total energ^ E = 1 in the system |IJ) with f3 satis- 
fying Q). There exists Sq > such that for any < e < Eq and for any path 
. . . <t_i<7o0"i ■ ■ ■ in the graph Z there exists a solution of (Qp and a sequence of times 
. . .t_it ti . . . such that the energies {5P of individual pendula satisfy 

\E*j fe) - !| < c Ve, and |£ CT (tj)| < Cy/e for a ^ a j: 

where C is independent of e. The times tj can be chosen so that 

< tj+i - tj < Ce- ir - 8 . (6) 

This theorem shows that, although the system (p} is near-integrable, so that for 
most (in the sense of Liouville measure) solutions the action stays close to its initial 
value for all time, there exist nevertheless solutions for which the action changes by 
O(l) no matter how small e is. In other words, the system exhibits Arnold diffusion. 
According to ((6]) the rate of this diffusion is polynomial. The bound in ((6]) is not 

2 In fact, any value in excess of the potential energy of an upside-down equilibrium works. 
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sharp, but it can be improved by a more careful tracing of the estimates in our 
example. In the general case, polynomial upper bounds for speed of diffusion for 
finitely-differentiable systems have been obtained in |Buj . 

The first example of Arnold diffusion was outlined in the well known paper of 
Arnold PQ. Bessi |Bs) (see also |BB| ) proved diffusion in Arnold's example by a vari- 
ational method, by considering the gradient flow of the Lagrangian action functional. 
John Mather |Maj used a somewhat similar approach to construct accelerating or- 
bits for time periodic mechanical systems on a 2-torus (see also |BT[IGTllDLS|lKa] ). 
References to the recent progress on Arnold diffusion goes beyond the scope of this 
paper and can be found in [KL2] . In the present paper we use a slightly different 
version of this approach, based on using the Maupertuis' principle. We construct 
the "diffusing" solutions as geodesies in a Jacobi metric, so that all these solutions 
have a fixed prescribed energy. These geodesies are constructed by concatenating 
geodesic segments which follow a prescribed itinerary. The construction is fairly 
similar to [KCTlKLl] . 

Anderson localization is an important example of energy (non)transfer (see 
[LTW] for a survey), still not very well understood. The role of Arnold diffusion 
for destruction of Anderson localization is discussed in |Baj . Probably the most 
popularized lattice model is the one introduced by Fermi-Pasta-Ulam in their sem- 
inal paper [FPU]. Although most small amplitude solutions in the FPU model do 
not exhibit energy transfer (see e.g. [HK] ) . proving the existence of solutions with 
energy transfer is an interesting open problem. Other physically significant lattice 
models are discussed in |FS W] . 

Understanding of the transfer of energy for Hamiltonian PDEs is one of emerging 
directions of research (see |Bo] ; a recent progress for the cubic defocusing nonlinear 
Schrodinger equation has been made in [I]). 




\ 
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>- : speed = 0(1) 



Advance from 123 to 234. 



Energy transfers from 1 to 2. 



Figure 4: One full step in the propagation of the "kink" . 
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1.1 A heuristic description of propagation 



In this section we give a purely heuristic picture of the physical motions showing 
Arnold diffusion. As mentioned earlier, we consider the periodic case Xi + <± = Xi as 
a representative example. 

Stage 1: transfer of energy. In this stage only three pendula: 1, 2 and 3 
governed by {1} with k = 1, 2 and 3 are "active" , while 4 "sleeps" upside-down 
(see Figured! left). The stage consists of many substages illustrated by Figures [SJ 
left. In each of these substages a small amount of energy is transferred from 1 to 
2. This transfer is somewhat similar to the one described in |KL1| for a metric on 
the 3-torus. Finally, in the last substage, 1 is left with just enough energy to climb 
upside down and to fall asleep there, while 2 rotates with speed O(l), as shown in 
the middle of Figure |H The same motion, viewed in the configuration space K 4 , is 
shown in Figure El 

The motion just described is similar to that in a slightly simpler example de- 
scribed in [KLT] [KL2] . 




Stage 2: advance. This stage is sketched in Figured] right and Figure[5l middle. 
At t = t\ three neighbors, say 1, 2, and 3, are in the bottom position. The middle 
pendulum 2 is running: ±2 = 0(1), while its two neighbors 1 and 3 have near- 
heteroclinic speeds close to the heteroclinic speeds y / —2V(x i ), i = 1,3 resp., at 
t = t\. The remaining pendulum 4 is up (see Figure |4] middle). As the time goes 
on, while 2 is spinning with speed 0(1), 1 rises to the top equilibrium, where it will 
sleep until further notice, while the sleeper 4 "wakes up" , i.e. falls from its perch, 
turning by 7r by the exact moment when 2 finishes a large integer number of full 
spins. By that moment, X3 makes a "gentle" turn by 27r, returning to the bottom 
position. In short, the accomplishment of this stage is the falling asleep of 1 and 
the awakening of 4. This is illustrated in Figured] right. We will call this stage 
the "advance" because of its similarity with the advancing caterpillar: a rear foot 
1 is placed on the ground, while the front foot 4 is lifted, ready to move. 
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The ending moment of the second stage is the beginning moment of the first 
stage described above modulo the shift of the index by 1. We have, in other words, a 
"traveling wave" - a (very) discrete analog of the kink in the sine-Gordon equation. 
However, in contrast to the standard traveling kink, ours can change the direction 
of its propagation arbitrarily, according to a prescribed itinerary. 

2 The proof of Theorem [D 

The full complexity of the problem is already seen in the case of 4 pendula, and 
we limit our consideration to that case. Now we restate the theorem in geometrical 
terms. The following is motivated by the heuristic outline of the energy transfer 
between the pendula: as mentioned before (see Figures |4] and [5]), we want the 
energy to pass from one pendulum (e.g. 1) to another (e.g. 2) in small increments 
over many steps, after which the pendula should change roles: the "giver" 1 "falls 
asleep" , while a "sleeper" (e.g. 4) "wakes up" thus enabling the "taker" 2 to become 
the next "giver". This is reflected in the following geometrical construction. First, 
we construct an itinerary (section 12. ip for the desired orbit. Then we reformulate 
the problem of existence such an orbit as a variational problem. Finally, we prove 
the existence of a solution to this problem in section 12.21 In the latter section we 
use two lemmas [3] and [2] stated (and proved) in sections [5] and [4] respectively. 

2.1 Constructing an itinerary 



string string 
q exchange ^exchange — advance exchange q exchange q 

^12,3 ^23,4 



Figure 6: An itinerary. 



To prove Theorem[T] given a path in the graph (see Figure[3]) we must produce an 
orbit that shadows this path. We describe the construction in the case of monotone 
energy transfer (i.e., to the right neighbor): (Tj+\ = o~j + 1 = j + 1 for all j G Z. 
Extension to the case of general path (. . . o§<j\<Ji . . . ) poses no new difficulties. We 
this consider an infinite sequence of codimension one sections in K 4 , grouped into 
finite strings, Figures [5] and [6] 

■■• (^12,3^12,3 ••• ^12,3) (^23,4) ^23,4) •••) ^23%) ••• ' CO 
^ ^ v. 

l->2 2^3 

where the sections and their spacings are defined according to the following rules: 

1. Section Y\ 2 3 , for example, is seen in FigureUl left. The subscripts 12 indicate 
that 1 and 2 exchange energy, and 3 is the "facilitator" : 

£12,3 :={&3 =0, x\ + x\ < e, \x 4 - tt| < yft). (8) 
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2. All sections within each string in ([7]) are integer translates of each other. For 
example, in the first string 1 — > 2: 

S 12 +3 = £13,3 + Mm k 12 , 3 , nf 2 , 3 , 1, 0) =: £f a>3 + nf 2 , 3 fc = 1, . . . , ^ 

To define the second string, we replace i H> i + 1 mod 4 in ((8J, setting 

S 23 , 4 := {^4 = 0, xl + xl< s, \ Xi --k\ (9) 

We then define £33 4 as a translate of £23,4 by integer multiples of 27r in each 
coordinate. Each section in the second string is an integer translate of the 
previous section by 

n 23 ,4 ! = 27r ( ' m 23,4: "23,4: l ) 

3. The neighboring strings are related via 

Center(j:l 3A ) = C 'enter (£^3) + 27r(i, m° 3 4 , 1, i) = Cenier(£^ 1 3 ) + r^ 3 4 . 

4. All translates are far apart: each 

\n k ..\>e~ 2r - 4 (10) 

5. The turns are gradual in the sense that the unit vectors e... — n k ./\n k .\ satisfy 

\e k+1 - e k ..\ < e 2r+ \ \e\. +1 - e N J \ < e 2r+i . (11) 
As mentioned above, treating a general itinerary poses no difficulties. 
2.2 A variational problem and its solution 

We note that the energy one solution of ([lj are the geodesies in the Jacobi metricd 



4 

1 - /XV{%i) - e/3(zi-i,Xi,x i+1 ,e)) ds, (12) 

i=l 

where x = (xi, x 2 , £3, £4), V(x) = — cosx — 1 < 0, ds is the Euclidean metric and 
indices of taken mod 4. 

With the sections having been defined in items 1-5 above, we now list the main 
steps of the proof of Theorem^ and fill in the details in the following sections. 

1. Defining geodesic segments. Let £o,£i be two consecutive sections in the 
chain of sections ([7]) and let Pi € Ex, i = 0,1. Centers of these sections 
differ by 2-7m with n being either (to, n, 1, 0) € Z 4 , or ( ^ , s, 1, |) with s G Z 
and satisfying (fTU|) . According to Lemma [3] from section [SJ there exists a 
connecting geodesic 7(po,Pi) of (TT2"]) which depends smoothly on its ends 
Po, Pi- At this stage the integer parameters either m and n or s are still free. 

3 Up to the factor v2, the square root is the speed of the energy one solution in the configuration 
space. 



dp(x) 



\ 
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2. Constructing a long shadowing geodesic. Consider a finite segment of N + 2 
sections from the sequence ([7]) . To simplify notation, we denote these sections 
by Sj, < i < N + 1 (N here is arbitrarily large). We also choose arbitrary 
points pi € Si . Later we will treat po , p jv+i a s fixed and pi , . . . , p^ as variable. 

According to the preceding item there exists a broken geodesic 7(^01 • ■ • iPn)> 
a concatenation of energy one orbits y(pi,Pi+i) of JT]). 

The length (in the Jacobi metric) of this broken geodesic 

L(pi, . . . ,p N ) = L(p Q ,pi) + . . . + L{p N -i,p N ). (13) 

is a function of break points pj\ we omit po and pn+i from the left-hand side 
since they will be considered as fixed. 

We will show that L has a minimum on the interior of its domain 

Si x ... x S Such an interior minimum corresponds to a true geodesic. We 
will thus establish the existence of a geodesic with a prescribed itinerary. 

3. The key step: existence of an interior minimum for A13\) . Let us consider 
two consecutive terms from the sum (|13l) : 

S(p) = L(pj-_i,p) + L(p,p j+1 ), p G S 2 (14) 

where Pj±i G Ej±i are fixed and p G Sj is variable. Without loss of generality 
we take Sj = S} 2 3 , given by To prove the existence of an interior 

minimum of (|13l) it suffices to show that the minimum of S(p) is achieved in 
the interior of Sj. To that end we first alter (3 in the Jacobi metric (TT21 by 
setting P — only in the cylinder x\ + x\ + x\ < e 2 which passes through 
the center of Sj (we do not alter (3 anywhere else). Before restoring /3 to its 
original form, we study the associated length 

S°(p) - L°( PQ ,p) + L°(p, P2 ), p G S, = S. (15) 

Once the properties of 5° are established (see (|16p and (|17p below), we will 
show that restoring /3 to its original form creates a minimum for S in the 
interior of S. Without the loss of generality, we take S = S} 2 3 as in (JH]), so 
that S a = S a (xi,X2,X4). 

We will first show that Sq = S°(xi,X2,X4 : ) is "nearly constant" in the first 
two variables and has a minimum near the "equator" X4 = n: 

|S°(xi,x 2 ,x 4 ) -5°(0,0, x£)\ <2ce r+2 - 5 , (16) 
S°(zi,Z2,7r±v^>S°(a:i,a;2,7r) + ! (17) 
for any (x\, X2, 7r) G S. 

4 For the future reference, we note that the triple may or may not be entirely in one string in 
the sequence J7J of sections. 
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Proof of (|16p and (|17p. By Lemma 1 from |KL1] we have, withp = (xi, x%, X3, 2:4): 



dL°(p_,p) _ dL a (p,p+) _. + 

d Xl Xi ' d Xl Xi ' 

where x _ (i) = (x^ , x^, x$, x^ ) is the energy one solution with the modified ft, 
connecting Pj-i to p, and where the differentiation is taken at the moment the 
solution passes through p. This solution exists by Lemma [3] below. We use a 
similar notation x + for the energy one solution connecting p with Pj+i- We thus 
conclude that 

— (p) = ±T _£+ i = 1,2,4; (18) 

this identitjH will allow us to analyze S 10 . Now, due to the fact that the perturbation 
P near p is removed the pendula are decoupled and the velocity is explicitly given 
in terms of energy distribution ([5]) 

\±i\ = y/2{Ei - V( Xi )), 

where Ei is the energy of the ith pendulum near p The estimate of (|18[) is now 
reduced to studying the difference of velocities. We have to consider two cases: in 
one, S, belong to the same string in (J7J (the case of "energy transfer"), 

and in the other, they do not ("the advance"). 



Case 1 — energy transfer. In this case all sections lie in the same string 
in |7]) - say, in 1 — > 2. The displacements 27m_ and 27m + are then of the form 
27r(m±, n±, 1, 0) with integer m± and n±. Assuming the integers to be positive (we 
can always assume them to be of the same sign) , we have 



±T > 0, x+ > 0, i= 1,2,3, 



(19) 



at the moment when E is crossed. 

The fact that the signs are the same for x + and x~ is of key importance because 
it provides a near-cancellation in (|18p for i = 1, 2, 3. Thus, we have 



dS° 

dxi 



^2(E- - V{x t )) - ^2{E+ - V( Xl )), i = 1,2. 



Note that if < A < B then y/B - y/A < \/B - A; this, used in $ZU§, gives 

dS° 



dxi 



< \ 2 E? -e; 



i = l,2. 



Now according to Lemma [2] and the assumptions (fTU|) and (| 1 1 1) we have 
\E+-Er\ <ce 2k+ \ 1 = 1,2,3,4 

and by (HJ) 

\d Xi S°\ <ce k + 2 , z = l,2. 



(20) 

(21) 

(22) 
(23) 



5 we do not differentiate by X3 since we only need to define S and S° on S C {23 = 0}. 
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Combining this with (j2"Tj) . we obtain (|16p . showing that S° is "flat" in the first two 
variables. 

To estimate dS° /dx4 note that (x^ , xj) and (x^ , x±) stay at the -^/i-neighborhood 
of the saddle (£4,2:4) = (tt, 0). Since the distance between sections is large (ITUl) and 
the phase velocity is bounded (by the choice of fixed energy), the duration of each 
stage is at least e~ 2k ~ 4 . This implies that (xj,xj) is at worst 0(exp(— e _1 ))-close 
to the unstable manifold 

y = U(x) = (x - tt) + 0((x - tt) 2 ), 

while (x~\,x\) is at worst 0(exp(— e _1 ))-close to the stable manifold y = —U(x) of 
the same saddle - all this at the moment when the solution crosses S. That is, 

xl = U{x) + Oiexpie- 1 )) (24) 

and 

xj = -U(x) + 0(exp(e- 1 )) (25) 

so that 

— = 2U(x) = 2(x 4 - tt) + 0((x 4 - tt) 2 ) + 0(exp(e- 1 )). (26) 
0x4 

Integration by x 4 gives (JTTJ). 



Case 2 — advance. In this case not all sections lie in the same string in (J7J; 
without the loss of generality, assume that Ej-i, £ are the last two in the string 
1—^2, while Sj+i is the first in the following string 2 — ► 3. In this case the 
displacement vectors are of the form ra_ = (m, n, 1,0), n+ = (i,s, 1, i). In this 
case we still have (|19|) . and following (f20|) and (|2T|) we obtain (fT6|) . Since the sign 
of £4 is unknown, we treat it separately, observing, as before, that (|24p and (|23|) 
hold. This implies ^ and thus (IT7)l . 

This completes the proof of ([TBI and p7|) in both cases. 

Proof of the interior minimium for S. Using the properties of S° and the 
positivity of /3 we now show that S has a minimum inside S. We do so for Case 1; 
the remaining case is treated almost verbatim. 

The boundary <9£ = d v Yl U c^X consists of the "vertical" and the "horizontal" 
parts (after possible reindexing of the coordinates): 

a ( ,S = {x\ +xl=e, \x4 - 7r| < Ve}, 
3ft E = {x 2 + x\ < e, |a;4 - ir\ = *Je}. 

A key observation we will use shortly is this: 

S(p) = S°(p) for all p G a„E. (28) 

Proof: We wish to show that the energy one solution r y(pj_x,p) with p lying on 
{x\ + x\ — e, X3 = 0} does not intersect the lens x\ + x\ + x\ < e 2 . To that end 
assume the contrary: the solution travels from one set to the other, taking some 
time t = t* > 0. Since the distance between the sets is > hyfe, while the speed 
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is < 2, the time of travel is t* > kV^- But ^3 > 1 f° r as l° n g as \ x a\ < 
Thus during time t*, a; 3 changes by the amount AX3 > \\fs, which means that the 
solution lies outside the lens x\ + x\ + x 2 < e, in contradiction with the definition 
oft*. This proves ([2"gjl. 

We now show that the restriction of S to each horizontal disk in E: 

D h ■= {x\ + x 2 2 <e, x 4 = 7T + h}, \h\ < y/e 

has a minimum in the interior of Dh- To that end, we first note a crucial fact that j3 
decreases S (as compared to S°) near the center Ch = (0, 0,n + h) of each Dh- Note 
that by the definition (Q} the infimum of /?(-, e) taken over the set x\+x^+x\ < e 2 /2 
is bounded from below by be r for some b > independent of e. Therefore, comparing 
the geodesic length in the original Jacobi metric with the truncated one, we obtain 

S{C h ) < S°(C h ) - einf /3(-,e) < S°(C) - be r+1 . (29) 

See [KL1] (proof of Lemma 4) for more details on this argument. 
On the other hand, by (fTS)) we have, for any p g dDh, \h\ < y/e: 

S a (C h )<S a (p)+ce r+2 - 5 . 

Combining this with (|29|) and (f28| we obtain 

S(C h ) < S(p) + ce r+2 - 5 - be r+1 < S(p), Vp e dD h . 

We showed that the minimum of S cannot be achieved on 9„£, and it remains to 
show that it cannot be achieved on c^E either. Estimate (fl7|) shows that S° has a 
pronounced minimum near the equator X4 = 7r. By the same estimate as we used 
for (|29|) . we have a two-sided result: \S°(x) — S(x)\ < e r+2 , which together with 
(ITTl) gives 

S(Xi, X2, 7T ± \fs) > S(xi,X2,7r). 

This proves that the minimum is achieved inside E. 

To complete the proof of the main theorem it remains to observe that the ex- 
istence of the the internal minimum for S implies the existence of the internal 
minimum for L(p±, . . . ,|>iv), as well as the existence of an internal minimum for an 
infinitely long sequence ([7]). The details can be found in [Lj. 

3 The pendulum Lemma. 

In this section we state and prove an auxiliary lemma which is used in the proofs 
of the main two lemmas in the following two sections. 

Lemma 1. For any \a\ < 1, \/3 — tt\ < 1, and T > there exists a unique solution 
x(t;T,a, ft) of x + sin a; = 0, satisfying x(0) = a, x(T) = f3 with the additional 
property 

a < x(t) < max{^,7r} for < t < T. (30) 

This solution depends smoothly on T, a and f3, and, moreover, as T increases from 
to co, the energy E — x 2 /2 + (—1 — cosx) decreases monotonically from 00 to 0. 
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A similar statement holds if either 



\a + n\< 1, |0| < 1, 



M < 1, |/3-2tt| < 1, 



or 



(3 - a > 2tt, 



with the solution confined ot the interval [— 1 — 7r, 1] in the first case, [— 1, 2tt + 1] in 
the second case and [a, [3] in the last case. In particular, in the last case f3 — a can 
be arbitrarily large. 



For any q,p e R 4 such that a = qi, f3 — pi satisfies conditions of Lemma\T\for each 
i = 1, . . . , 4, then for any T > there is a unique solution X{t, T, g, p) of the system 
i31\) satisfying X(0) = q and X(T) = p. Moreover, there is a unique T > such 
that energy of this solution is one. 

We first prove the Corollary. 

Proof of the Corollary. Lemma[T] applies to each of the n equations in (|3"TT) . by the 

assumptions of the Corollary. That is, for any T > there exists a unique solution 
Xi(t) — Xi(t; T, qi,Pi) of Xi+smxi = with the desired boundary conditions. Now as 
T increases from to oo, the energy Ei(T) of each solution decreases monotonically 
from oo to 0. There is therefore a unique T = T* with £™ =1 i^(T*) = 1. In addition, 
T* is a smooth function of q,p, by an application of the implicit function theorem. 



Proof of the Lemma. We concentrate on the first case; the remaining ones are 
essentially the same. We write the equation of the pendulum as the system 



We are seeking a solution starting on the line x = a and ending at t = T on the 
line x — j3. The condition j3 > ir — 1 imposes a lower bound on the initial velocity. 
This leads us to consider the ray OM of initial data on the line x = a, see Figure 
[3 where O is the point whose solution crosses the x-axis at x = n — 1. 

Consider now the ray OM carried by the flow of (|32|) for time T, where T is 
fixed arbitrarily. Consider the set It on the ray OM whose points enter the strip 
S = {tt — l<x<7r + l}at some time t <T and do not leave it before t = T. We 
claim: For any T > 0, the set It is an interval, and its image ip 7 (It) under the 

6 Since we chose to concentrate on n = 4 pendula, we formulate the lemma for this case, although 
the proof carries over verbatim for an arbitrary n. 



Corollary 1. Consider the uncoupled systen^ 



Xi+smxi=0, i = l,...,4. 



(31) 



Q.E.D. 




(32) 
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Figure 7: The slope of the image of the vertical interval is positive. 



flow at t = T is a curve with a positive slope, connecting the two boundaries of the 
strip S, Figure^ This claim implies the statement of the Lemma. 

If T is small, i.e. if the solutions are fast and thus lie above the separatrix, the 
result is obvious. However, for T large some solutions starting on It "turn around" , 
as in Figure [7] (B), and the proof requires a little care. 



Let zq = (a,yo) € It', this, we recall, means that (i) z T = tp T zq 6 S, and (ii) 
once z* enters S it does not leave S before t = T. Consider the linearization of (|32p : 

l'\ ., (33) 
rj = — (cosx) £, 

where x is the solution of (|3"2j) with the chosen initial condition. Note that the 
solution Q — (£,r)) of (13"3"|) with £(0) = (0, 770)1 ryo > is a tangent vector to the 
image curve ip T (OM) at the point <p T zq. It suffices, therefore, to prove that C(^) 
lies in the first quadrant. This equation is changing type from elliptic to hyperbolic, 
and thus is it not a priori clear that the solution may not execute an unwanted 
rotation, violating positivity of the slope rj/t;. The idea is to show that during the 
"dangerous" elliptic stage while x(t) < tt — 1, r]/^ remains positive. 

To that end, let r 6 (0, T] be the time of entrance of the solution y'zo = 
(x(t),y(t)) into S, so that 

x(t) € [a, ir - 1] for t € [0,t] 

and 

x(t) e [tt- 1,tt + 1] for t e [t,T] 



We will first show that > at t = t. 

The key idea is to to observe that the vector z = (x, x) rotates clockwise faster 
than the vector £ = (£,77); this will be shown shortly. Since the slope of z is positive 
at t = r, the same will then be true of £(r) as desired. We claim: 

^Mrf(zA <0 , wnenev er » = 2. (34) 
dt\xj dt\tj 1 x £ y ' 
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To prove (|34[) we carry out the differentiations and use the equality of the slopes to 
reduce the inequality to an equivalent one: 



or, using and ([55]): 



x 2 < ^ <U ' 



-x since — y 2 — cosx£ 2 — rj 2 

x~ 2 < e < 



Using the equality of slopes, this reduces to 

sin a; 



> cosx, 
x 

which holds true thanks to x £ (0,7r). We conclude: since y(r)/x(r) > 0, and 
V (r)/ar) > 0. 

We show that the slope of £ remains positive for the remaining time [r, T] . 
During this time \x — ir\ < 1 and thus cosx < 0. Hence the linearized vector field 
crosses into the first quadrant, and since £(r) lies in that quadrant, it is still 
there at t = T. 

The monotonicity of E = E(T) also follows from (i) the positivity of the slope of 
the image curve ip T (It), and (ii) from the fact that the curve moves "to the right" as 
T increases. Indeed, if we increase T, the point of intersection of ip 1 \It) with x = f3 
moves down, to the curve with the same x — /? but with the smaller y, i.e. with 
the smaller value of E = y 2 + (— 1 — cos (3). Furthermore, if T is small, the velocity 
must be large: (3 — a — x(T) — x(0) < (max y) T, so that max y > — > oo as 
T — > 0. On the other hand, if T is sufficiently large, then all solutions starting on 
It must pass arbitrarily close to the saddle (tt, 0) and thus must have energy close 
to that of the saddle, i.e. to -E sa ddic = — (— 1 — cos7r) = 0. 

The proof of Lemma [T] is complete. 

4 The Hyperbolic Lemma. 

By Corollary Q] of Lemma [I] given any points q, q' € R 4 with their coordinates 
Xi and x\ lying [—1, 1] mod 2-k or in [it — 1,7T+ 1] mod 27r, there exists a solution 
X(t; q, q', T) of (f5T|) with e = which travels from q to q' in time T. By the same 
corollary, there exists a unique T(q,q') for which the total energy of the solution 
X{t;q,q',T(q,q')) is one: 

4 X 2 
^Ei=l, where E i = ^- + V(Xi). (35) 

i=l 

We thus associate with the energy one solution (of (HJ with e — 0) connecting q and 
q', the energy vector 

E(q,q') = f (E 1 ,E 2 ,E 3 ,E 4 ); 
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according to the Lemma [I] this vector is uniquely determined by the endpoints q, q' . 



Lemma 2. If two pairs of points: qi,q[ and qi,q'i in K 4 satisfy the condition^ 

\q'k-qk\>e- 2r -\ k = l,2, (36) 
and i 

|e 2 -ei| <e 2r+5 , where efc = , g f ~ gfc fc = l,2, (37) 

m - ik\ 

then the energy vector E of the connecting solution of with (3 = satisfies 

|E(g 2 ,4)-E( gi ,g0l<£ 2r+5 - (38) 
Moreover, there exists a constant C such that for all q, q' with \q' — q\ > 1 we /icwe 



^X(0:q,q',T(q,q')) 



< C; (39) 



/iere i/ie notation ^ is used to emphasize that the q-dependence enters X in two 
places - one through the boundary condition, and the other through T(q,q'). 

Proof. Statement (|38|) follows from the proof of Lemma [1] the main difficulty is in 
proving (j3"9]ffl To prove (|39|).we expand its left hand side: 

£-X{0; q, q\ T(q, q')) = d g X(0; q, q', T) + d T X(0; q, q', T) ■ d q T{q, q'), (40) 

where T = T(q,q') is to be substituted after the differentiations on the right-hand 
side. We will now estimate each of the summands on the right-hand side separately. 

Estimate of d q X(Q\ q, q',T). The proof of Lemma [T] shows that each component 
Xi of X depends on the boundary conditions xi , x\ and T onlj@, but not on xj , x'j 
with j i. This implies that the matrix d q X(0;q,q',T) is diagonal, with the 
diagonal entries dXi(xi, x\, T)/dxi. But this derivative is simply the slope of the 
image of the line x = x' under the map f~ T , where ip 1 is the phase flow of the 
pendulum equation. The argument of Lemma [1] shows that, because of the shear 
in the phase velocity field, this slope is always bounded once T exceeds a fixed 
constant. It remains to prove the upper bound for the last summand in (|40[) . 



Estimate of 9yX(0; g, q',T) ■ d q T(q,q'). We will first show that this term is ex- 
pressible via the first factor alone, thus reducing the number of estimates needed. 
Note that dxX ■ d q T is a square matrix, the product of the row matrix &tX with 
the column gradient matrix d q T{q,q'). To prove the lemma, it remains to show 
that each entry 

|d r X i (0;s ij iJ J T).a Xi r( 9) g / )|<C J (41) 

for \q' - q\ > 1. 



7 Since we chose to concentrate on n = 4 pendula, we formulate the lemma for this case, although 
the proof carries over verbatim for an arbitrary n. 

8 this estimate can be strenghened: C can be replaced by C/\q' — q\, but we do not need this in 
our proof. 

9 We recall the notation q = (x\ , 12, £3, £4). 
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Some identities. Let 



/ * 1-, X , I , and K = s Tk s . (42) 

J Xi {2{E i -V{x)))m) ^ 

We will show that 

5 Xi T(g, g / ) = ^" 1 M(0;a ;i) ^ ) T) ) (43) 



thus reducing (|41[) to an equivalent inequality 

IA- 1 M(0; Xi,^,T) drXj{Q;xj,x' j: T)\ < C. (44) 

Heuristically, one expects that |9T^i(0;a;i,x-,r)| < cT^ 1 . Indeed, Xi(0] x i7 x' i7 T) 
is the y-coordinate of the intersection, in the (X, X) -plane of the line {X = Xi] 
and the curve tp~ T {X = x'^, where (p t is the phase flow of the pendulum equation. 
Now because of the shear in the phase flow, one expects the line It — <p~ T {X = x^} 
to form angle at most cT -1 with the trajectories. Thus the point It H {X — Xi] is 
expected to move with speed < cT -1 , suggesting that indeed |<9t^(0; x^x'^T)] < 
cT^ 1 . We carry out a precise proof by an alternative, purely analytical method 
(which ultimately reduces to the same estimates) . Namely, we will use the following 
identity: 

which, together with (|43l) . is proven in a separate section below. 



Estimate of K 1 . Since Y,j =1 Ei = 1, we have j < Ei < 1 for some i, and thus 
for some C we have 

js—l = f 5 dx [< dx = 

1 J Xi (2(E i -V(x)))^- J Xi ^2{E t -V{x)) 

so that 

K-^C^Kj) <K^<CT. (46) 

Estimate of <9t^(0; X*, x' { ,T). We consider two separate cases: (i) |x^ — Xi\ < 2ir 
and (ii) \x' { — Xi\ > 2ir. 

1. In case (i), an estimate of (|45|) is easier done geometrically, as follows. Consider 
the graph y — Ut{x) of the time T-preimage of the line x — x' in the phase plane 
of the pendulum. Let y — U (x) be the graph of the unstable manifold of the saddle 
(it, 0); by a standard hyperbolic argument, the flow takes the line exponentially close 
to the unstable manifold: \Ut{x) — U(x)\ < e~ cT , for \x\ < 7T, and, moreover, the 
motion of the line becomes exponentially slow: 

\±U T (x)\<e- cT : for |x| < tt; (47) 
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here T is greater than a fixed positive constant because of the assumption \q'— q\ > 1. 
But U T (xi) = X i (0;q,q',T) and (gTJ) gives 

\d T X t (Q;q,q',T)\ < e~ cT , for \x[ - x t \ < 2tt. (48) 

This completes the proof of (|44p . and thus of the Lemma, in case (i). 

2. In case (ii) we have x\ — xi — 2irn,i + r, < r < 2tt with integer n ^ 0. We will 

use (|45|) to prove (j44)l , to which end we need an upper bound on Ki . From (|42"T) we 

have 

i da; da; n, , . 

*f 1 > rn / > / 77"— TT^H = *ET, (49) 



where c~2n/ + w 2 /2. Now the number of revolutions rii > T/TE t where T(Ei) 
is the time of one full revolution: as 

T Ei = P dx < s/2 r d \ = ^(Ml+Vl + ED-lnEt). 

J -* y2(E. t + 2 sin 2 x/2) J ° V E i + ^/^ 

Substituting this into (j4"9)) we get 



v ^ Ei EfTEi . Ei(ln{l + ^/T+Ei)- In Ei) 
Ki<c—< c^r- < ci . 

Finally, we substitute this estimate into (1431) : 



IdTX^x^x^T)] < c — = < - 



Together with ([4T)f this proves (|44p . The proof of the lemma is thus complete. 



Proof of the identities (|43p and (|45p. The energy of the solution X(t; Xi, x^, T) 
is a smooth function of Let Xi,x' i} T; we denote this energy by E(xi, a;-, T), so that 

X(0-Xi,x'i,T) = ^{Eixux'^-Vixi)). 
Differentiating by T we get 

dl y /2(E(xi,x' i ,T) - V(Xij) 

To estimate the numerator, we differentiate the identity 

T=p - dx (51) 

J., v /2(75(a; l ,a;^T)-y(a;)) 

with respect to T and solve for dE/dT, obtaining 

we^-v^ Y 1 <52) 
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Substituting this into (|5U|) proves (|45l) . To prove the remaining identity (|4"3")l . we 
recall that T(x, x') is the time which gives energy one to the solution: 



^2E(x k ,x' k ,T(x k ,x' k )) = 1. 



(53) 



fc=i 



Differentiating this by gives 
^(ii,a:J,r) 



<9x,- 



T=T(x,x') 



dT(x,x') j^ dE(x k ,x' k ,T) _ Q 



<9x,; 



fc=l 



(54) 



The above sum, according to (|52|) . can be replaced by — St=i ^ft = — -^j solving 
for dT/dxi gives 



dT(x,x') _ dE(xi,x' i ,T) 



dxi 



dxi 



T=T(x,x') 

To estimate the last derivative, we differentiate the identitj0 (|5"Tj) by ^ 







1 



dE(x t ,x' % ,T) 



y/2(Ei - V( Xi )) 4 (2(Ei-V(x))f/ 2 d Xi 



k: 



dEixux'^T) 



K, 



dx z V /2(S, - V( Xi )) 

Substituting this into (1551) results in the proof of (|4"3| : 

3T(x,x') ! g5j 

-K 



dx t 



^2{E, - V(xi)) 
The proof of the two identities is now complete. 



(55) 



5 The connection Lemma. 

The following lemma is the building block in the construction of shadowing geodesies. 

Lemma 3 (Existence of geodesic segments ). There exists eo > such that for 
all < e < So the following holds. Consider any two sections, which we denote by 
So and Si, in the itinerary ([?}), such that the vector 2i:n connecting the centers of 
S , ^i satisfies \n\ > -. This vector is of the form 2it(m ) n, 1, 0) or 27r(i,s, l,i) 
( with integer m,n or s ), modulo possible translation in the index. Then for all po € 
So, pi G Si there exists a geodesic j{po,Pi) in the Jacobi metric H2\) connecting 
Po with p\ , and depending smoothly on po and p\ . 
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Figure 8: Towards proof of Lemma [31 



Proof. 

1. We first define a section So shown in Figure as follows. By Lemma [1] 
there exists a (unique) solution X (t) of (f3"T|) of energy one, connecting c = 
Center(T,o) with c\ — Center{Yi\). We now define the section Sq as the 
codimension one disk in R 4 at the distance £3 from Co and perpendicular to 
the initial direction eo = -Xo(0)/|Xo(0)|: 

So = {q : (q- c ) • e = |g - (c +£3e )| < £*}, 

where • denotes the usual dot product. Analogously, we define the section Si 
near ci, except for reversing the sign of e\. 

2. For any pair of points go G So, 9i G Si we consider three geodesic segments 
(in the metric (fT2|) ): y(po>?o)) 7(<Zo>9i) an d 7(gi,Pi), along with the velocities 
of the associated solutions of ([T]) u^, ur, iol, wr as shown in Figure leff 2 ! . The 
lemma will be proven once we show that there exists a pair goi Qi , smoothly 
dependent on po, p\, for which 

Vl = vr and w L = wr. (56) 

To that end we first list the properties of each of the three geodesic segments. 

3. Since radius of Eo is e 1 / 2 , we have \qo — j?o I — O(e^), which is small compared 
to the injectivity radius of the metric (fl2j) . By standard arguments from 
differential geometry (using the smoothness of solutions of the ODEs and the 
implicit function theorem) we conclude that 

VL=v y^ — —\ + r OL (po,q ,e), \vol\c 1 = 0(e*); (57) 
\qo - Pol 

here ^o = 3y/2 is the speed at Co (by the energy restriction (|3"5l and the fact 
that the center of the lens corresponds to three pendula at the bottom, and 

10 Here and below E t = E(x i ,x' i ,T). 

11 Since we chose to concentrate on n = 4 pendula, we formulate the lemma for this case, although 
the proof carries over verbatim for an arbitrary n. 

12 It should be noted that vl, Vr, wl,wr all depend on qo, q\. 
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one at the top.). A similar estimate holds for the right end: 

wr = v j^ — +r 1R {q ll p l ,e), \r XR \ c ^ = 0(e*). (58) 

\Pi - 9i| 

4. The intermediate segment 7(^0, 9i) avoids the lenses, and thus Lemma ap- 
plies; in particular, the C 1 -bound (|39]) holds, implying that 

VR = X (to)+roR(q ,qi,e), \r 0R \ c i <Ce^. (59) 

and 

w L = X (h) + ri L {qi,pi,e), \r 1L \ c i < Ce^ . (60) 
Here (i = 0, 1) is the time when Xo(t) intersects the section Si. 

5. We will prove the existence of the pair qg , qi satisfying (|56l) by applying the 
implicit function theorem. To that end let v denote the orthogonal projection 
of v £ R 4 onto R 3 D S; wc will also treat q a G R 4 as an element of S C K 3 , 
denoting it by qo € R 3 . To prove (|55|) it suffices to prove that the projected 
equations 

vl = v R and w L = w R (61) 

hold, since if (|6Tj) hold, then the remaining components orthogonal to S must 
match as well by the conservation of energy. Substituting the estimates (|57jl , 
([55]). (HH and (HOI) into J5B]) we obtain, after projecting onto R 3 D 5: 

«o, 1 = r (po,qo,qi,£), v - , = n(q , qi,pi,e) (62) 

w-po\ \pi-qi\ 

as the equivalent matching conditions, with the C^-small remainders: 

|fi| c i < Cei 

In arriving at (p32"j) . we made use of the fact that X (t ) — 0(si), as follows 

from the choice of Sq to be orthogonal to X o (0) (so that X n (0) — 0) and the 
the fact that to = 0(e3). 

6. To apply the implicit function theorem, instead of the variables gi we introduce 

Qq = v - Qi=v - 



\ql-po\ \pi-m 

Qi £ R 3 . Expressing 

9o = Po + v \qo -po\Qo, 9i = Pi - Vo\pi - qi\Qi 
and substituting into (|61|) . we obtain 

Qo = Ro(po,Qo,Qi,£), Qi = Ri(Qo,Qi,pie)- 

Introducing Q = (Qo,Qi) € R 6 and R = (R ,Ri) we rewrite the matching 
condition (|56p in the final form 

Q = R{Q,Po,Pi,e), 
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where < Ce^ . It is important to observe that R is defined (at least) 

on the entire ball \Q\ < |, independent of e. Thus for all sufficiently small s 
there exists a unique solution Q depending differentiably on the parameters 
Po, Pi- 

This completes the proof of Lemma [3] 
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